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Neutrino Oscillation Search at MiniBooNE 

Z. Djurcic* for the MiniBooNE Collaboration a t 
a Columbia University, New York, NY 10027, USA 

This article reports the status of a ~ * u e oscillation search in MiniBooNE (Booster Neutrino Experiment) 
experiment. If an appearance signal is observed, it will imply Physics Beyond the Standard Model such as the 
existence of light sterile neutrino. 
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1. INTRODUCTION 

In recent years, the solar-neutrino pQ, reactor- 
neutrino [2], atmospheric- neutrino [3], and 
accelerator-neutrino [4] experiments have con- 
firmed the existence of neutrino oscillations. 
These results implied the existence of two inde- 
pendent Am 2 regions, with Am 2 ~ 8 x 10 _5 eU 2 



in the solar, and with Ar 



3 x 10" 3 eU 2 in 



the atmospheric sector. The Standard Model in- 
corporates two separate oscillation regions with 
three known neutrino flavors: v e , v^, and v T . 
However, an unconfirmed evidence for neutrino 
oscillations came from the LSND [5] experiment 
with Am 2 at ~ leU 2 value. The discovery of 
nonzero neutrino masses through the neutrino os- 
cillations has raised a number of very interesting 
questions about the neutrinos and their connec- 
tions to other areas of physics and astrophysics. 
One question is whether there are sterile neutri- 
nos that do not participate in the standard weak 
interactions. This question is primarily being 
addressed by the MiniBooNE experiment. The 
MiniBooNE experiment will confirm or refute the 
LSND result with higher statistics and different 
sources of systematic error. If the LSND neutrino 
oscillation evidence is confirmed, it will, together 
with solar, reactor, atmospheric and accelerator 
oscillation data, imply Physics Beyond the Stan- 
dard Model such as the existence of light sterile 
neutrino [6]. While LSND observed an excess of 
v e events in a beam, MiniBooNE is a — ► v e 
search. 
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2. THE EXPERIMENT 

MiniBooNE is a fixed target experiment cur- 
rently taking data at Fermi National Accelera- 
tor Laboratory. The neutrino beam is produced 
from 8.89 GeV/c protons impinging on a 71 cm 
long and 1 cm diameter beryllium target. The 
target is located inside a magnetic focusing horn 
that increases the neutrino flux at the detector 
by a factor of ~5, and can operate in both neg- 
ative and positive polarities for v and v running. 
MiniBooNE collected approximately 6 x 10 20 pro- 
tons on target (POT) in neutrino mode. This 
data sample is currently used in the neutrino os- 
cillation analysis. Mesons produced in the target 
decay-in-flight in a 50 m long decay pipe. The 
neutrino beam is composed of from K + /ir + — > 
/i + + decays. The neutrino beam propagates 
through a 450 m of a dirt absorber before entering 
the detector. There is a small contamination from 
v e . The processes that contribute to the intrinsic 
u e in the beam are /i + — > e + v e v^, K + — > 7r°e + ^ e , 
and K° L -> ir ± e ± v e . Early in 2006, MiniBooNE 
switched the polarity of the horn to select nega- 
tive sign mesons. 

The MiniBooNE detector is a 12.2 m diame- 
ter sphere filled with 800 tons of pure mineral oil. 
The center of the MiniBooNE neutrino detector is 
positioned L = 541m from the front of the beryl- 
lium target. The vessel consists of two optically 
isolated regions divided by a support structure 
located at 5.5 m radius. The inner volume is 
the neutrino target region, while the outer vol- 
ume forms the veto region. There are 1280 8- 
inch photo-tubes (PMTs) pointed inward provid- 
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ing 10% coverage, and 292 outward-faced PMTs 
in the veto region. Data analysis require a fidu- 
cial volume cut at 500 cm from the center of the 
tank to ensure good event reconstruction, result- 
ing in a 445 ton target region. The outer volume 
serves as a veto shield for identifying particles 
both entering and leaving the detector. The rate 
of neutrino candidates per proton delivered to the 
target was constant (1.089 x 10~ 15 v/POT) over 
the period of data acquisition, as demonstrated 
in Figure [T] Protons delivered to the MiniBooNE 
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Figure 1. Calibrated neutrino rate per proton on 
target. 



trino induced events are identified by requiring 
that the event is observed during the beam spill, 
have fewer than 6 veto PMT hits, and greater 
than 200 tank PMT hits. These simple cuts yield 
a cosmic ray rejection that is better than 1000:1. 

The MiniBooNE detector is calibrated using a 
laser system, cosmic rays and neutrino interac- 
tions. The calibration procedures cover the full 
energy range from 50 to 1000 MeV using differ- 
ent event types. Stopping cosmic ray muons are 
used to calibrate the energy scale for muon-type 
events and measure the position and angle recon- 
struction resolution, when their path length can 
be identified. This is accomplished with a scintil- 
lator hodoscope on the top of the detector, com- 
bined with scintillation cubes at various positions 
within the detector volume. Analysis of the cos- 
mic muons indicates that the energy and angle 
resolution of muons is 10% and 6.4°, respectively. 
Observed Michel electrons from muon decay are 
used to calibrate the energy scale of electron-type 
events at the 52.8 MeV Michel endpoint, with 
12% resolution. Reconstructed n° events provide 
another electron-like calibration source. The pho- 
tons that are emitted in tt° decays span a consid- 
erable range to over 1000 M eV. The 7r° mass de- 
rived from reconstructed energies and directions 
of two 7-rays has a peak at 136.3±0.8 MeV. The 
resolution on the reconstructed 7r° mass peak is 
20 MeV. This is in an excellent agreement with 
the 135.0 MeV expectation, providing a check on 
the energy scale and the reconstruction over the 
full energy range of interest for the v e appearance 
analysis. 



target are monitored by two toroids located in the 
beam line. The main MiniBooNE trigger is an 
accelerator signal indicating a beam spill. Every 
beam trigger opens a 19.2 fis window in which 
all events are recorded. Other triggers include 
a supernova trigger, a random strobe trigger for 
beam-off measurements, a laser calibration trig- 
ger, cosmic muon triggers, and a trigger to record 
neutrino events from the NuMI beam line. 

The interaction point, event time, energies, and 
the particle tracks are recorded from the times 
and charges of the PMT hits in the detector. Neu- 



3. NEUTRINO FLUX, CROSS SEC- 
TION, AND DETECTOR MODELING 

The flux modeling uses a Geant4-based simula- 
tion of beam line geometry. Hadron production in 
the target is based on Sanford-Wang parametriza- 
tion of p — Be cross-section, with parameters de- 
termined by a global fit to p — Be particle produc- 
tion data. Simulated neutrino flux has an energy 
distribution with a peak at ~ 0.7GeV. There- 
fore, the average LjE v ratio is ~ 0.8km/GeV 
compared to LSND's L/E v ~ lkm/GeV. 

The NUANCEv3 [7] event-generator simulates 
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interactions in the detector. The cross section 
model describes the various neutrino interaction 
processes on CH 2 , which include the Llewellyn- 
Smith free nucleon quasi-elastic cross section, 
the Rein-Seghal resonant and coherent pion pro- 
duction cross section, a Smith-Moniz Fermi gas 
model, and final state interactions based on n- 
Carbon scattering data. 

The MiniBooNE detector is modeled using an 
extended GEANT3-based simulation. An ex- 
tended light propagation model of the detector 
describes the emission of optical and near-UV 
photons via Cerenkov radiation and scintillation. 
Each photon is individually tracked, undergoing 
scattering, fluorescence, and reflection, until it 
is absorbed. The response of the electronics to 
the photoelectrons resulting from photons hitting 
PMT's is simulated, with the final output be- 
ing digitized waveforms simulating the charge and 
time channels of the electronics in a form identical 
to that used for data. 

4. THE BACKGROUNDS IN THE AP- 
PEARANCE SIGNAL 

The signature of an oscillation event is the 
v e + n — > e + p reaction. The backgrounds in 
the oscillation analysis are divided in two main 
categories: intrinsic v e events in the beam, and 
ZAj events that are mis-identified as v e events. 

The beam that arrives at the detector is almost 
pure with a small (0.6%) contamination of v e 
coming from muon and kaon decays in the de- 
cay pipe. The v e from /i-decay are directly tied 
to the observed interactions. Taking into ac- 
count a small solid angle subtended by the Mini- 
BooNE detector, the pion energy distribution can 
be determined from the energy of the observed 
events. The pion energy spectrum is then used to 
predict the v e from /^-decay. A second source of 
v e originates from if e 3-decay. This component 
is constrained using high energy charged cur- 
rent quasi-elastic (CCQE) events, that originate 
primarily from kaon decays. 

Mis-identified events are 7r°, A-decays, and 
CCQE events. Most 7r° are identified by the 
reconstruction of two Cerenkov rings produced by 
two decay 7-rays, as shown in Figure^ However, 




Figure 2. (Preliminary) Mass distribution of ir° 
candidates obtained using a two ring likelihood 
fit. The data are the black points, the black line 
is Monte Carlo simulated neutrino interactions, 
while the red line is the subset of these events 
which have at least one pi in the final state. 



the decay of 7r° can appear much like primary 
electron emerging from a v e charged current in- 
teraction if one of the gammas from the decay 
overlaps the other, or is too low in energy to be 
detected. Over 99% of the NC tt° are rejected in 
the appearance analysis. In addition to its pri- 
mary decay A — > tt N, the A resonance has a 
branching fraction of 0.56% to the 7 N final state. 
The 7-ray may mimic an electron from v e interac- 
tion. The rate of A production in neutral current 
interactions can be estimated from the data, us- 
ing the sample of reconstructed tt° decays. Most 
Vy, events can be easily identified by their pene- 
tration into the veto region when exiting muons 
fire the veto, or by muons stopping in the inner 
detector and producing a Michel electron after a 
few microseconds. 
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5. SIGNAL SEPARATION FROM THE 
BACKGROUND 

Particle identification (PID) is performed by 
different algorithms that use the difference in 
characteristics of the Cerenkov and scintillation 
light associated with electrons, muons, protons or 
i s, as shown in Figure [3l These algorithms in- 
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electrons following the primary interaction, and 
requiring the veto to have less than six hits to en- 
sure there is no cosmic muon contamination and 
the tank to have greater than 200 hits to sup- 
press decay electrons from cosmic muons. The 
events are then fit under the single electron and 
muon ring hypotheses and a ratio is formed with 
the resulting likelihoods. Events with likelihood 
ratios favoring the muon hypothesis are rejected. 
The event is then fit with two ring fits, both with 
the mass free and fixed to the nominal tt° mass. 
In the case from Figure 02 separation between ir° 
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Figure 3. Cartoon showing how particle identifi- 
cation acquire an input from the Cerenkov rings. 
From the top to the bottom, the rings are: an 
ideal particle, an electron, a stopping muon, and 

7T°. 



elude a maximum likelihood method and boosted 
decision trees [8] . Figure [4] shows the logarithm 
of the likelihood ratio formed from fitting neu- 
tral current ir° candidates under a single electron 
ring hypothesis and a two ring hypothesis where 
the kinematics are fixed to give the nominal ir 
mass. Electron-like events should have larger pos- 
itive values, while 7r°'s should have negative val- 
ues. The events are selected by requiring no decay 



Figure 4. (Preliminary) The logarithm of the like- 
lihood ratio formed from fitting neutral current 
7T° candidates under a single electron ring hypoth- 
esis and a two ring hypothesis. The data are the 
black points, the black line is Monte Carlo sim- 
ulated neutrino interactions, while the red line is 
the subset of these events which have at least one 
7T in the final state. 



and v e candidates is achieved. Similar PID sep- 
aration is performed between electrons and other 
types of detected events. An example of a boosted 
decision tree where the muon/electron separation 
was measured with cosmic ray muons and associ- 
ated electrons is given in [5]. The PID removes 
- 99.9% of v^, CCQE interactions, - 99% ef- 
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fective 7r° producing interactions, and preserves 
a high efficiency for z/ e interactions. An opti- 
mized PID is expected to allow a small contami- 
nation of an oscillation signal with mis-identified 
tt° (~ 83%), A-decays (~ 7%) , and CCQE 
events (~ 10%). 



6. BLIND ANALYSIS STRATEGY AND 
CROSS-CHECKS 

MiniBooNE is conducting a blind analysis in 
order to complete an unbiased oscillation search. 
That means that the region where the oscillation 
z/ e candidates are expected is closed for the anal- 
ysis. In the example given in Figure 21 the mass 
from the fit is required to be greater than 50 
MeV/c 2 . The likelihood ratio of the single elec- 
tron fit and the fixed mass two ring fit is required 
to favor the 7r° hypothesis. The latter two re- 
quirements are for blindness, in order to keep the 
region with potential oscillation v e candidates out 
of the analysis reach until the final set of selection 
cuts is formed and the PID is optimized. In prac- 
tice, the MiniBooNE analysis cannot use its own 
data to verify PID algorithms with an oscillation- 
like v e data sample. However, an important cross- 
check of electron event reconstruction and par- 
ticle identification comes from NuMI events ob- 
served in the MiniBooNE detector. Neutrinos 
produced by the decay of mesons moving along 
the NuMI beamline and in the vicinity of the 
NuMI target may reach the MiniBooNE detec- 
tor. NuMI events consist of v e , u^, n^, ir° and 
A over the range of energies relevant to the ap- 
pearance analysis. A kaon and non-kaon fractions 
in NuMI Monte Carlo are extracted from a fit to 
the data. Such Monte Carlo is then compared 
to our Monte Carlo prediction. The NuMI sam- 
ple therefore allows an independent check of the 
MiniBooNE PID algorithms performance in iso- 
lating v e . Other important checks are performed 
with either MiniBooNE or external data. A com- 
plete list of cross-checks is given in Table [T] 



Table 1 

The MiniBooNE analysis is verified by different 
experimental cross-checks for each event class rel- 
evant to v e appearance search. 
K+ HARP [H3], External Data 

MiniBooNE Data 
K° E910 [EE], External Data 

MiniBooNE Data 
fi MiniBooNE Data 

7T° NuMI, MiniBooNE Data 

Other (A, etc) NuMI, MiniBooNE Data 



7. REMAINING STEPS IN THE OSCIL- 
LATION SEARCH 

When the final set of the analysis cuts is deter- 
mined and associated systcmatics evaluated, the 
data sample that potentially contains the oscilla- 
tion candidates will be un-blinded. The compo- 
sition of the final sample will be predicted by the 
MiniBooNE Monte Carlo simulation, with Monte 
Carlo sample filtered through the same set of PID 
cuts. If MiniBooNE confirms the LSND result, an 
excess of events in the data distributions when 
compared to Monte Carlo will be observed. The 
final event sample will be evaluated using a x 2 
function 

x 2 = ^(o, - mc^)- 1 ^ - Pj), (i) 

»=1 

applied to the data and Monte Carlo energy dis- 
tribution of the oscillation candidates. Oi is the 
number of observed events in an energy bin i. Pi 
is the Monte Carlo prediction that takes into ac- 
count oscillation parameters (sin 2 29, Am 2 ). The 
covariance matrix Cy account for the statistical 
and systematic uncertainties in the energy bins. 
Systematic errors are associated with neutrino 
flux, neutrino cross sections, and the detector 
model. 

The flux prediction has the uncertainties cor- 
responding to the production of tt, K, and Kl 
particles in the MiniBooNE target. These un- 
certainties are quantified by a fit to external data 
sets from previous experiments on meson produc- 
tion. The cross section uncertainties are eval- 
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uated by continuously varying underlying cross 
section model parameters in the Monte Carlo con- 
strained by MiniBooNE data. Uncertainties on 
the parameters modeling the optical properties of 
the oil in the MiniBooNE detector are constrained 
by a fit to the calibration sample of Michel elec- 
trons. The uncertainties are currently being eval- 
uated. 

8. CONCLUSION 

The MiniBooNE experiment will confirm or re- 
fute the LSND oscillation signal with approxi- 
mately 6 x 10 20 protons on target collected for 
the analysis. The oscillation appearance analysis 
is underway, with current work on the systematic 
error evaluation that combines the errors from the 
v flux, v cross-sections, and the detector model- 
ing. 
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